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Abstract 

We use the density matrix renormalization group (DMRG) to map out 
the ground state of a XY-spin chain coupled to dispersionless phonons of 
frequency u>. We confirm the existence of a critical spin-phonon coupling 
a c tx u for the onset of the spin gap bearing the signature of a Kosterlitz- 
Thouless transition. We also observe a classical-quantum crossover when the 
spin-Peierls gap A is of order oj. In the classical regime, A > u, the mean-field 
parameters are strongly renormalized by non-adiabatic corrections. This is 
the first application of the DMRG to phonons. 
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The spin-Peierls state [|l| in quasi-one-dimensional materials has attracted renewed at- 
tention lately because of its discovery in the organic series |J (BCPTTF) 2 X , with 
(X = AsFq, PFq), and in the cuprate compound [|J CuGeO^. Our own interest in this field 
stems from earlier findings on non-adiabatic effects or quantum lattice fluctuations on the 
transition temperature T c of the closely related (through the Wigner- Jordan transformation 
HI) electron-phonon problem. These generally tend to decrease T c || f| and the order pa- 
rameter || and can even destroy it for spinless fermions when the phonon frequency gets 
appreciably larger than the gap |§. The existence of this adiabatic-quantum crossover has 
thus far remained only a theoretical concept. Moreover, there are predictions of a Kosterlitz- 
Thouless (K-T) transition in the quantum regime || which yet need to be confirmed. With 
the advent of the density matrix renormalization group (DMRG) |H| , numerical calculations 
of the ground state properties of large scale systems has now become a reality. Although 
phonons have not yet been approached by this method, we will show the feasibility of such 
calculations. 

We propose to study a XY-spin chain whose magnetic interaction depends on the bond 
length. The reason is twofold. First, the uniform spin chain enjoys an exact solution. It 
can be mapped into a half-filled non-interacting spinless fermions (pseudo-electrons) chain 
of width 2 J using the Wigner- Jordan transformation Q. Second, the model contains the 
essential elements for a spin-Peierls transition, that is coupling to inter-molecular motion. 
We have chosen a dispersionless vibrational spectra for the molecular motion along the chain 
direction. This is quite acceptable as it is well known that it is the phonons at the zone edge 
that couple to the spins. This also circumvents the contribution of the hydrodynamic modes 
to the molecular displacements which increase logarithmically with chain length. This is 
especially annoying if one wishes to study the thermodynamic limit of the system as one 
has to go to long chain lengths. We start from the following Hamiltonian 



H = Y, + [ J + a (&} + b e - b\ +1 - b t+1 )] (SfSf +1 + SjSj +1 ) } 



in which u is the vibrational energy (or frequency 
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annihilation (creation) operator for a vibration on molecule £, J is the magnetic interaction, 
a = gj\f%j where g = (dJi : £ + i/d£) is the standard spin-phonon interaction in the classical 
paper of Bray et al. [|lj, and is the local XY spin of value \. 

The expected crossover and the nature of the different regimes warrants further ex- 
amination at this time. In the two-cutoff analysis ||, the natural energy cutoff for the 
pseudo-electrons is the bandwidth for characteristic energies (e.g. the gap A) larger than u. 
This is the classical regime in which the adiabatic approximation prevails. But for character- 
istic energies less than u, there is a new cutoff uj and the pseudo-electron-phonon coupling 
generates instantaneous backward-scattering g\ w — 8ct 2 /(u;7rvF) and Umklapp g^ ~ —g± 



couplings (in the renormalization group notation [TlJ JT2[ ||13||), where vp is the Fermi ve- 
locity of the pseudo-electrons. This is the quantum regime in which the interactions are 
unretarded and quantum fluctuations impregnate the ground state. It turns out that, for 
spinless fermions, the local character of the backward scattering cancels out. It is only the 
non-local character (wave number dependence) that can give rise to non-trivial effects. The 
renormalization equations || show that unless the non-local contribution to #3 has the 
right sign and has a bare (initial) value larger than a certain threshold, the Umklapp pro- 
cesses are irrelevant (g^ renormalizes to zero) and the quantum system is gapless. If the 
threshold condition is met, the Umklapp processes and the vertex function grow to infinity. 
This signals the onset of a gap in the excitations. The transition is believed to be of the K-T 
type [|l5j 0. We thus propose to look for a threshold value a c for the appearance of a gap 
and for the signature of the Kosterlitz-Thouless transition. The formula derived by Baxter 
T6fl for the order parameter P is of the form P oc A -1 exp (— A -1 ) where is A oc (T c — T) 1 ^ 2 



for a thermodynamic transition. For a ground state transition, identifying T with g\ oc a 2 
suggests A oc (a 2 — a 2 ) 1 ^ 2 . We shall refer to this last form as the effective K-T coupling. 

We now turn to the use of the DMRG. Its inherent difficulty in dealing with vibrations 
is the infinite dimension of the Hilbert space for phonons. The most direct way around this 
problem is to truncate the space. This is obviously unsatisfactory whenever the average 
molecular displacements involve too many virtual phonons rii = \(be)\ 2 . A crude estimate 



for this can be obtained by equating the increase in elastic energy to half the decrease in 
electronic energy U. For single electrons, as in the polaron problem, one gets is ri£ w U/2u 
which obviously precludes any serious study of the deep adiabatic region uj <^U [HJ with a 
limited phonon space. In the case of a collective state, such as the spin-Peierls (or Peierls) 
modulated state, the same analysis for non-interacting pseudo-electrons leads, in a mean-field 
type analysis, to the criterion ng ph A 2 ln( J/ A)/2umJ. One gains a factor A ln( J/ A) / irJ over 
the single-electron situation, allowing a deeper incursion into the adiabatic region uj <C A. 
In our case, we typically kept ri£ less than 3 and u was as small as A/25 under the best 
conditions. 

In the DMRG procedure a set of internal sites are coupled, via the inter-site part of the 
Hamiltonian, to environmental blocks. The DMRG proposes an iteration algorithm for the 
growth of the environmental blocks and thus, for the size of the system. Our adaptation 
of the method, for each set of parameters (a/J,u/J), was the following. At each iteration, 
we considered a single site coupled to the environment blocks generated by the previous 
iteration. This central site had two spin states and M v vibration states, equal to anywhere 
from three (at small gaps) to twenty (for the larger values of a). M v could be estimated by 
requiring that the local vibration subspace properly describe a coherent state of n v phonons, 
that is with an error much less than the allowed DMRG truncation error. The total number 
of sites in the chain was chosen to be even such that, with open boundary conditions, the 
chain would relax to a unique broken-symmetry state. Consequently, the number of sites of 
the bordering environment blocks differed by one, there being one long block and a one-site 
shorter block. We kept M& state in each block, a number varying from 60 to 120 (for the 
smaller gaps), and the corresponding matrix elements of all operators coupling the central 
site to the bordering sites of each block. We targeted the four lowest energy states in the 



superblock using the Davidson-Liu [TJj algorithm, two in each of the total spin subspaces 
having S z = 0,1. This way, we were able to reliably access the spin gap energy on the ansatz 
that, in the spin-Peierls state, the ground state is a singlet S = state and the spin gap is to 
an S = 1 state. These four target states were then used to project out the density matrix pij 



of the two new blocks formed by the direct product of the spaces of the left (right) block and 
the added site. The z component of the spin being a good quantum number, we used it to 
fragment the block spaces into more manageable sizes. The density matrix was diagonalized 
and the M& highest energy eigenstates were kept. The truncation error Tr p — J2i=iPu, 
which was used as a minimal error indicator, was kept at a few times 10~ 5 or less. The 
block sizes thus increased by one and the chain length, by two sites at each iteration. The 
total length varied from 100 to 600 sites at the smaller gaps. At the longer chain lengths, 
the numerical accuracy on the gap greatly suffered and finite size scaling arguments were 
used to analyze the data. The error on the gap is estimated to be 10% under these strained 
conditions. 

Fig. 1 shows typical values of the gap plotted in such a fashion as to bring out the 
characteristic K-T form A oc (a 2 — a^) _a5 exp (— b(a 2 — a;?) -0 ' 5 ). The signature is clear and 
there is indeed a critical spin-phonon (pseudo-electron-phonon) coupling a c for the onset or 
long-range order. Its values are obtained by least-squares fitting and are plotted in Fig. 2. 
This phase diagram identifies, for the first time, the existence of a power-law relationship 
a c ~ uj ' 7 . The average molecular displacement u = 
law although the fit is less reliable as u varies much more slowly. Fig. 3 hints to the 
probable existence of a power relationship between A and u which one should expect if both 
parameters follow the K-T formula. We could not define one unambiguously. 

Fig. 3 points to the identification of a crossover between two regimes. One, on the left of 
Fig. 3, is characterized by a rapid increase of A/au and the other one, on the right, shows 
a rather flat dependence A/au ~ 1. As a matter of comparison, the mean-field solution [|J 
predicts A/au = 2. It is thus tempting to identify this latter region, which has a mean- 
field like behavior, as the classical or adiabatic regime described earlier. Fig. 4 confirms 
this assignment. It shows a crossover from a classical region, with a gap of the mean- field 
form A = A Q exp(— AJu/a 2 ), into another regime, the K-T one of Fig. 1. Moreover, this 
crossover occurs at the expected value A « uj (at the position of the arrows in Fig. 4). 
This crossover is also shown in the phase diagram (Fig. 2). The classical region is obviously 
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also seems to obey a similar 



limited to A < uj. We have found the gap to saturate at a value very near J for large uj. This 
is why the adiabatic regime is restricted to uj < J in Fig. 4. The gap prefactor A Q and the 
leading coefficient A in the exponential are strongly frequency dependent as shown in Fig. 5. 
We find that A extrapolates to 0.44 and A Q to 1.3 J at zero frequency (calculations cannot be 
performed at uj — > since the number of virtual phonons goes to infinity). These values are 
close to the expected mean-field ones of 7r/8 and 4J/e, respectively |§. The difference is 
significant, however, in view of a least-squares error of order 2% for these parameters. It is 
perhaps our linear extrapolation from finite frequencies which is improper. This conclusion 
is borne out by the plot in Fig. 5 of the plateau value of A/cm in Fig. 3 (an average 
value in the classical region). It clearly shows that the classical regime does not satisfy the 
mean- field condition |l| A/au = 2 and that the extrapolation to uj — > is surely far from 
linear. It is however likely that this ratio can reach the value of 2 and that our results are 
consistent with the mean-field limit at uj = 0. The parameters A D and A in the classical 
regime expression for the gap are seen in Fig. 5 to be strongly renormalized by non-adiabatic 
phonon contributions. Moreover, there seems to be a crossover, within the classical regime, 
between two behaviors. For uj < .05J, the renormalization of A G and A are identical. This 
suggests a common origin to the renormalization, i.e. that it is the magnetic interaction J 
which changes at these smaller frequencies. For J > uj > .05J, these two parameters evolve 
differently, the extra renormalization of A coming seemingly from vertex corrections to a. 
The reason for this peculiar behavior is unclear. 

We were also able to monitor the change in phonon frequency in the classical regime. 



We find, as Wu et al. JT9[ did, that the renormalized phonon frequency uj r decreases as 
A = 8a 2 /iTUjJ does. For instance, we find at the crossover uj/A = 1 that uj t /uj « 0.5, 0.7 for 
uj = .025, 0.3 or A « 0.3, 1.3. This renormalization is stronger than the one Wu et al. found 
for spin 1/2 fermions. 

In conclusion, we have not only verified the predictions of a Kosterlitz-Thouless transition 
and of a classical-quantum crossover in the spinless-fermion-phonon problem but we have 
also been able to produce a quantitative phase diagram of the spin-Peierls system and deduce 
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a power law dependence of the critical spin-phonon coupling on the phonon frequency . We 
have also shown that the DMRG can be used advantageously on phonon problems be they 
only dispersionless in this paper. 
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FIGURES 



FIG. 1. Semi-log plot of the spin gap A as a function of the reciprocal of the effective 
K-T coupling, for two values of the phonon frequency. The lines, as well as the quoted values 
of a c , are least squares fits to the formula of Baxter in the quantum region. 

FIG. 2. Phase diagram of the system. The full line is a fit of the Kosterlitz-Thouless 
transition line. It has the power-law dependence a~ 2 oc oo~ lA . The dashed line straddles the 
classical-quantum crossover Za ~ lu and ends at uj ~ J. 

FIG. 3. Ratio of the spin gap A to the lattice modulation amplitude u times the spin- 
phonon coupling a as a function of the square of the effective K-T coupling, for two values 
of the phonon frequency. The arrows point to the axis each data set refers to. 

FIG. 4. Semi-log plot of the spin gap A as a function of the reciprocal square of the 
spin-phonon coupling a, for two values of the phonon frequency. The straight lines are least 
squares fits to the mean-field like expression in the classical regime. The arrows point to the 
axis each data set refers to and also indicate where the classical-quantum crossover occurs. 

FIG. 5. Evolution of the parameters A D and A of the mean-field type expression for the 
spin gap A as a function of frequency u, normalized to their zero frequency values A o (0) 
and A(0) which were linearly extrapolated from the data at finite frequencies. Also shown 
is the frequency dependence of the average value of A (u)/au in the classical region, where 
a the spin-phonon coupling and u is the average molecular displacement. This latter ratio 
should equal 2 in a mean-field solution. The arrows point to the axis each data set refers to. 
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